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Phases in optical lattices vs. Coulomb frustrated HTc cuprates 
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Fermionic atoms in 2D optical lattices and electrons in HTc cuprates may both be described 
by the Hubbard model. However, if Coulomb frustration is responsible for the striped phases 
in 2D cuprates the phase diagrams will differ markedly. Two representative scenarios are de- 
scribed by a simple stripe model without phase separation and a mean field model with phase 
separation in the absence of Coulomb frustration. When Coulomb frustrated both models display 
antiferromagnetism (AF) and stripe phases with d-wave superfluidity, whereas neutral atoms in 
optical lattices will only do so in the stripe model. Radii and densities of the various phases in 
harmonically confined optical lattices are calculated for the two models and have very different Mott 
plateaus and density discontinuities. Observation of antiferromagnetic, stripe and superfluid phases 
in density and momentum distributions and correlations from time-of-flight experiments is discussed. 
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I. INTRODUCTION 

Ultra-cold atomic Fermi-gases present a new opportu- 
nity to study strongly correlated quantum many-particle 
systems. Optical lattices provide a periodical lattice po- 
tential in which the atoms are described by the tight- 
binding approximation, which enables access to the Hub- 
bard model with interactions tuned by Feshbach reso- 
nances. Most importantly we gain direct observation 
into many controversal issues in strongly correlated high 
temperature (HTc) superconductors in a controllable way 
where interactions, densities, temperatures, etc., can be 
tuned. Recent experiments with optical lattices have 
measured momentum distributions and correlations, and 
have found superfluid phases of Bose and Fermi atoms 
Pi 0, 0]> Mott insulators pj, and band insulators [a, @|. 

HTc cuprates, however, suffer long range Coulomb re- 
pulsion between doped electrons or holes (if localized) 
whereas neutral atoms in optical lattices are not likewise 
frustrated. Coulomb frustration inhibits phase separa- 
tion and may also be responsible for stripe formation 
in cuprates which again affects d-wave superconductiv- 
ity (dSC). By studying the phases in optical lattices we 
can directly observe whether phase separation, stripes 
and dSC occur in the Hubbard model without Coulomb 
frustration for general density, temperature, interaction 
strengths, etc. 

Many models have been applied for investigating HTc 
such as the 2D Hubbard and t-J Hamiltonians in dif- 
ferent approximate versions with very different results 
(see [24( and references therein). Extending with next- 
nearest neighbor hopping and interactions or long-range 
Coulomb interactions further complicates the problem. 
Early Hubbard mean field (MF) calculations predicted 
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phase separation between an antiferromagnetic (AF) 
phase with density n — 1 and a paramagnetic phase for 
U^7t or a ferromagnetic phase for U^7t [7|. Phase sepa- 
ration (PS) in the Hubbard model and the closely related 
t-J model is still controversial despite intense investiga- 
tions. In the t-J model PS is found for large J/t but for 
the more realistic smaller values of J/t the field is divided. 
In the Hubbard model without next-nearest neighbor in- 
teractions {t' = 0) exact diagonalization [§| and Monte 
Carlo results [§] find PS whereas dynamical mean field 
[Io| and variational cluster perturbation theory do 
not. In dynamical cluster approximations PS is found 
for finite t' [H]. 

A MF ground state with stripes was found by Za- 
anen and Gunnarson These stripes have a hole 

density of one per site whereas density matrix renor- 
malization group calculations find a stripe density 
of 1/2 in accordance with experiments fl5|. However. 
Green's function Monte Carlo calculations p(| find only 
weak signs of stripes. Various cellular dynamical MF [17j 
and cluster calculations [18| find competing antiferromag- 
netic and dSC phases. The long range Coulomb forces in 
cuprates suppress phase separation, possibly leading to 
mixed phases with short range density fluctuations such 
as stripes [20]. The atoms in optical lattices are electri- 
cally neutral and the two phases of different densities are 
therefore not forced to mix by Coulomb forces to e.g. a 
stripe phase but can remain in separate bulk phases. 

The many different results indicate that the ground 
state is delicately balanced between nearly degenerate 
phases. Slight approximations or changes in hopping and 
interaction parameters can change the ground state phase 
dramatically. Thus, two decades after the discovery of 
HTc, the ground state and the mechanism of pairing re- 
mains undetermined in cuprates as well as Hubbard and 
t-J models. 

The purpose of this work is rather to study qualitative 
differences in confined optical lattices from simple cal- 
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culations for the 2D Hubbard model in which Coulomb 
frustration play a central or no role for inhibiting phase 
separation, stripe formation and HTc. To this purpose 
we specifically investigate in section II two simple rep- 
resentative models with and without PS, their different 
density distributions and phases. We investigate the ef- 
fects of Coulomb frustration in section III and show that 
it is difficult to distinguish the two models in HTc ma- 
terials because phase separation is inhibited and stripe 
phases appear in both cases. Likewise dSC appears in the 
stripe phases as discussed in section IV. Yet in confined 
optical lattices as described in section V, the absence of 
Coulomb frustration leads to very different phases, den- 
sity and momentum distributions and correlations for the 
two models. 



II. EQUATION OF STATES FOR THE 2D 
HUBBARD MODEL 



In the tight binding approximation spin 1/2 fermions 
are described by the Hubbard model [21j 



{ij)a i 



(1) 



where di CT is the usual Fermi creation operator (a =|, J,), 
n i(J = a^Ofo the density, and (ij) denotes nearest neigh- 
bors. Only 2D square lattices at zero temperature will be 
studied here. U is the on-site usually repulsive interac- 
tion and t the nearest-neighbor hopping parameter. The 
model can be extended to include next-nearest-neighbor 
hopping (f') 0, [22j and longer-range hopping and inter- 
actions. We will mainly discuss the t' = case where 
the phase diagram is symmetric around half filling. It 
is then sufficient to discuss only the doped case where 
x = l-n>0. 

As mentioned in the introduction, various approxima- 
tions and numerical results to the Hubbard model lead 
to a variety of phases as ground state solutions. Long- 
range Coulomb repulsion present for the case of elec- 
trically charged electrons/holes further complicates the 
phase structure and will be included in a later section. 

In the following two subsections two simple represen- 
tative models, a MF and a stripe model, are investigated 
with and without PS respectively. 



A. Phase Separation in Mean Field models 

MF theory provides a first impression of the phases 
competing for the ground state and has the advantage 
that it is computationally simple as compared to more 
complicated theories. The MF equations for the Hub- 
bard model are standard and we refer to, e.g., Refs. 

The energy densities can be calculated within the 
Hartree-Fock approximation for the paramagnetic (PM) , 
ferromagnetic (FM), antiferromagnetic (AF) and other 



phases. In Fig. 1 we show an illustrative example for the 
energy density vs. filling fraction for U = 6t and in the 
symmetric case t' = 0. 

At low density n <C 1 the ground state is that of a 
dilute paramagnetic (PM) gas with energy 
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(2) 



In units where the lattice spacing is unity (a = 1) this 
energy per site is also the energy density, and the density 
is the site filling fraction. 

Near half filling e PM = -(4/7r) 2 t + U/4, which be- 
comes positive when the repulsive interaction exceeds 
U/t > 64/7T 2 ~ 6.5. The PM phase is then no longer 
the ground state. In a state with only one spin the an- 
tisymmetry of the wavefunction automatically removes 
double occupancy and the repulsive term Un 2 /4 in Eq. 
([2]) disappears. Such a ferromagnetic (FM) state always 
has negative energy £fm < for n < 1 and is a candi- 
date for the ground state. AF, linear AF [23] and stripe 
phases are other competing candidates. 

For the 2D Hubbard model it is, however, necessary to 
extend MF to mixed phase scenarios. When U^7t the 
ground state of the MF Hubbard model undergoes tran- 
sitions from an AF at half filling to a mixed AF+PM 
phase for |x| up to a finite (U dependent) value where- 
after a pure PM phase takes over. [7] For larger U the 
phase diagram is more complicated with a pure as well 
as mixed FM phases between the AF and PM phases. A 
finite next neighbor hopping term t' makes the phase di- 
agram asymmetric around x = 0, extends the AF phase 
and changes the phase diagram considerably. 

At half filling n — 1 the ground state is an AF. Near 
half filling < x -C 1 the MF equations and AF energy 
can be expanded as 



£ A F = ~J[1 + -^X 
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(3) 



The concave dependence on x signals phase separation 
into a mixed phase of AF and PM by the Maxwell con- 
struction as shown in Fig. 1. The PS extends from the 
AF phase with density n = 1 to a PM phase with U de- 
pendent density n s — 1 — x s ~ 0.6. The double tangent 
determines the PS energy eps within MF. 

The difference in energy between the PS and the 
PM/FM phases sets the energy and temperature scale 
for the AF correlations, which is related to the pseu- 
dogap, which in turn is related to the spin gap as [25| 
T* ~ A s /2. We expect that the spin gap is of order the 
energy difference between the AF phase and spin uncor- 
rected PM or FM phase, whichever has lowest energy. 
Thus T* ~ (e PS - e PM )/2 - (1 - x/x s )J/2. The AF 
spin correlations have a long but finite length scale so 
that the Mermin- Wagner theorem is not violated. 
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Figure 1: 2D Hubbard MF energy densities in the case U = 6t 
are shown for the PM, AF and the Maxwell construction for 
the phase separation (PS). Also a stripe energy is shown (see 
text). 

B. Stripe model 

Stripes have been discovered in low-energy magnetic 
neutron scattering in doped cuprates at incommensu- 
rate longitudinal and horizontal charge and spin wave 
numbers [l5j, e.g. Q c = (27r/a)(0, ±x/x s ) and Q s = 
(7r/a)(l, 1 ± x/x s ) for the longitudinal charge and spin 
wave numbers respectively. Here, the hole filling in the 
stripes is half filled x s = 1/2 for low dopings x < 1/8 
but at larger doping it increases linearly to filled stripes 
x s = 1 such that x/x s — 1/4 remains constant. There- 
fore the charge (spin) density stripes appear with peri- 
odic distance d (2d) depending on doping as d — ax s /x. 
The stripe distance decreases with increasing doping un- 
til x > 1/8, whereafter the stripes remain at a distance 
d — 4a. The stripes act as anti-phase domain walls and 
the spin density wave therefore has periodicity twice the 
length of the charge-density wave. 

Diagonal, horizontal/vertical, checkerboard stripe so- 
lutions have been found in a number of models. In 
MF models Zaanen and Gunnarson [13J found stripes of 
hole density x s = 1 that are vertical or horizontal for 
U /t^2> — 4 and diagonal otherwise. In DMRG calcula- 
tions [3| stripes with hole density x s = 1/2 are found in 
agreement with experiments flB |. 

Insight into pairing and stripe formation can been ob- 
tained by adding or removing one, two or more parti- 
cles/holes in the half-filled AF. When the hole moves 
a distance of I lattice distances the kinetic energy is re- 
duced by t/l 2 but the AF is frustrated by an energy ~ I J. 
Solving such a simple linear string-like model for the one- 
hole problem in a 2D AF for J^t leads to a hole localized 
within a distance I ~ (i/J) 1 / 3 with energy 

£l = -At + (3J 2/3 t^ 3 . (4) 

Within the WKB approximation one obtains /3 = 



(37r/4) 2 / 3 ~ 1.8. Self consistent Born approximations 
[H find similar hole energy e x /t = -3.28 + 2.16( J/t) 2 ' 3 
valid when J^QAt. 

Because two nearby holes of opposite spin can move 
through the AF ordered lattice without upsetting the AF 
order they are not localized. If the pair has an energy 
that is lower than that of two separate holes e 2 > 2si—e , 
the odd-even difference can as in nuclear physics be taken 
as an effective pairing energy A = (erj + £2)/2 — £i. 

Adding more holes one might expect that more pairs 
form, eventually creating a molecular Bose-Einstein con- 
densate (mBEC) as in the BCS to BEC crossover in 3D 
ultracold atomic traps [2fJ, [27]]- On the square lattice 
such pairs have d-wave symmetry and mBEC will be a 
dSC [2a |. However, a number of numerical calculations 
as well as experiments find that the holes form periodic 
stripes rather than forming a mBEC. The equi-distance 
of the stripes indicate that they "repel" each other even 
in calculations without Coulomb frustration. As a con- 
sequence the chemical potential must decrease with hole 
density or equivalently e" = (d 2 e/dx 2 ) x= o > 0, and we 
can expand the stripe energy density as 

e= - J-e!x+^e"x 2 + 0{x 3 ) 7 (5) 

for small doping. It is the second derivate e" of the energy 
density at small doping that distinguishes the phases. It 
is positive for the PM, FM and stripe phases, zero for PS 
but negative for the unstable AF. 

The stripe phase is a specific ordered mixed AF and 
PM phase and is a continuous transition between the two 
pure phases as function of density. Similar mixed phase 
solutions are believed to occur in neutron star crusts be- 
tween nuclear matter and a neutron gas f29ll . and possibly 
also between quark and nuclear matter [301 . In both cases 
Coulomb energies add complexity to the mixed phases by 
ordering them into structured crystallic phases. 



III. COULOMB FRUSTRATION 

Long range Coulomb interactions prevent phase sepa- 
ration into two bulk phases of different charge density. 
They also inhibit the formation of localized holes, pairs 
and stripes. The phase diagrams of cuprates with such 
Coulomb frustration and optical lattices without may 
therefore be very different. 

Coulomb frustration in cuprates has been discussed in 
connection with stripes (see e.g. [25l. l3ll|). Generally the 
CDW hole pairs and stripes are energetically less favor- 
able due to long range Coulomb repulsion. In the follow- 
ing we shall consider the stripes as rods with charge less 
than the surrounding phase and calculate the additional 
Coulomb energy of such structures. 

Coulomb energies have been calculated for structures 
of various dimensionality D and volume filling fraction / 
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Here, the charge density difference between the two 
phases is Ap = ex s /a 3 for the stripes and the volume 
filling fraction is / = x/x s . The filling fraction is also 
the inverse of the stripe distance in lattice units. The 
dimensionality is D = 3 for spherical droplets or bub- 
bles, D — 2 for rods and tubes and D = 1 for plate- 
like structures. The diameter of the spheres, rods or the 
thickness of the plates is of order the distance between 
layers 2R ~ a ~ 4. The stripes are rods in a 2D plane 
but are embedded in a 3D layered structure. For D = 2 
the expression in the square bracket of Eq. J6]) reduces 
to [ln(l//) — 1 + /]. The logarithm originates from the 
Coulomb integral J dz/z along the rod length z, which 
is cutoff by other rods at a length scale I ~ a\/J. For the 
cuprates we furthermore reduce the Coulomb field by a 
dielectric constant of order e ~ 5. The resulting Coulomb 
energy of stripe or rod-like structures D = 2 is 



Ec ^^/[ln(l//)-l 



/] 



(7) 



Energy costs associated with the interface structures 
are usually added. Such surface energies are difficult 
to calculate for the stripes because their extent is only 
a single lattice constant. In principle they are already 
included in the stripe models. We will therefore just 
add the Coulomb energies given above. However, the 
Coulomb energies and the energy of the systems as a 
whole, may be reduced by screening and hole hopping 
into the AF whereby R increases but Ap and / are re- 
duced. 

Inserting numbers e 2 /tlc — 1/137, e = 5, x s — 1/2, 
a = 4 we find that e c ~ 150 meVx[ln(l//) - 1 + /]. 
In comparison the energy gain by changing phase from 
an AF to a stripe or PM phase increases with doping 
as (saf — £ s /pm) ~ Jx. The linear coefficient is J = 
(ex — 3 J/2) for the stripe phase but considerably smaller 
for the PM phase J ~ 50— 100 meV according to Fig. 1. 
The Coulomb energy of Eq. |[7j) thus dominates at small 
doping due to the logarithmic singularity and therefore 
the AF phase of density n^,l is preferred. 

The AF phase is the ground state as long as the 
Coulomb energy of Eq. (0 exceeds Jx corresponding 
to doping less than 



xaf — x s exp 



8 Jae 
7r x s e 2 



(8) 



Inserting the above numbers and J ~ J we find xaf — 
0.1 which is within range of the observed \xaf \ — 0.03 for 
hole doped and xaf — 0.15 for particle doped cuprates. 
In MF the particle-hole asymmetry arises from the next- 
nearest neighbor hopping t' ~ — 0.3i and leads to an AF 
phase extending from half filling up to a particle doped 



density n > 1. [3] Note that it is not taken into account 
that with increasing doping the AF density approaches 
that of the PM, and eventually the charge difference Ap 
between the AF and PM and the resulting Coulomb en- 
ergy become sufficiently small that stripes are favored. 

In the above picture the incommensurate stripe phases 
at small doping arise due to Coulomb frustration when 
t' = 0. At larger doping the stripes approach and will 
eventually affect each other. Experimentally the stripes 
undergo a transition from an incommensurate to a com- 
mensurate phase at x ~ 1/8 corresponding to a stripe 
periodicity of four lattice spacings. 



IV. SUPERFLUIDITY 

Since phase separation and stripe phases are sensitive 
to Coulomb frustration, we can ask the related question: 
does superfluidity occur in 2D optical lattices with Fermi 
atom as in HTc cuprates or is Coulomb frustration re- 
quired? The answer depends on the mechanism behind 
HTc which twenty years after its discovery is still not 
well understood. Like the ground state phases discussed 
above, the various models disagree about the origin of 
dSC and its dependence on the Hubbard parameters, 
doping and the influence of stripes and Coulomb frus- 
tration. In certain models [25j], stripes are a prerequisite 
for dSC. In the following a simple semi-phenomenological 
model will be discussed based on similar ideas that pair- 
ing take place on stripes in an AF background. 

Standard BCS superconductivity is s-wave and re- 
quires an attractive pairing interaction. HTc supercon- 
ductivity is d-wave and is believed to arise in the cor- 
related state of the Hubbard model with purely repul- 
sive forces. The effective pairing must be a result of 
highly collective effects since the two-body interaction 
is strongly repulsive which inhibits onsite s-wave pairing 
but favors next-neighbor d-wave pairing between oppo- 
site spins due to super-exchange. 

An instructive example where superfluidity arises is 
the ID Hubbard model in a staggered or AF background 
magnetic field. At low magnetic field it is a normal 
metal as the standard ID Hubbard model [32] whereas 
for strong field it becomes superfluid for x < (2 — y2) 
[33T ]. The staggered magnetic field leads to a linearly in- 
creasing string potential between holes of opposite spin 
and constitutes a strong effective pairing potential. How- 
ever, the stripes observed in experiments have antiphase 
domain walls. The staggered magnetic field is therefore 
of opposite sign on the two sides of the stripe and no net 
staggering magnetic field results. No superconductivity 
is thus expected on the ID stripes if detached from the 
surrounding 2D AF phase. 

Instead, effective pairing correlations in two-leg lad- 
der models are found transversely to the stripes [25| . 
Two features of superfluidity are important: that the 
electrons or atoms must pair and they must condense. 
These two features limit the critical temperature for weak 
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and strong pairing respectively. In the weak limit the 
critical temperature is exponentially suppressed as in 
standard BCS T c ~ exp(-l/g) (for 2D see Eq. p3]l). 
where g is the product of the pairing attraction and the 
level density. In the strong limit it is argued [25[ that 
T c = H 2 n s (T = 0)/4m* where n s is the superfiuid den- 
sity and m* the effective mass of the particle. The super- 
fluid density is limited by the density of holes and thus 
vanishes with decreasing doping. 

A simple model will now be constructed that incor- 
porates these ideas for the inhomogeneous stripe phase 
in a semi-phenomenological way. The weakly interacting 
limit, where the effective pairing and T c vanish, is ap- 
proached when the AF correlations dissolve and the PM 
phase appears. To estimate the pairing in this BCS limit 
we assume for simplicity that the pairing take place be- 
tween nearest neighbor holes and, since it is driven by 
the AF correlations, that its effective strength V is pro- 
portional to the volume of the AF. Specifically, we shall 
assume that it vanishes when there is only one AF stripe 
between each half filled stripe, i.e. 



V = V (l-2x/x,). 



(9) 



Here, the effective pairing strength must exceed Vq > 2t 
and Vo^7.35t in order that holes can form s-wave and d- 
wave two-body bound states respectively [lH]. In contrast 
the onset of superconductivity occurs for Vq > 2t and 
Vq > respectively. 

In the strongly interacting limit at low hole density 
we ensure that the superfiuid density vanishes with hole 
density by explicitly multiplying the density of states by 
the hole density. 

Generally, we assume the pairing takes place between 
pairs in the stripes, although pairs can tunnel between 
stripes. The derivation of the gap equation for holes then 
follows that of standard BCS on a lattice (see e.g. f22|). 
The pairing is expected to be d-wave with a gap A q = 
2A[cos(<7 :c ) — cos ( (?,,)l since s-wave is suppressed by on- 
site repulsion [22j, [3J|. The resulting dSC gap equation 
for the critical temperature T c is 



1 

V 



x 1 v^r / \ , ^2 tanh(e q /2T c ) 

^ [cos(q x ) - COs(q y )} ^3 . (10) 



x, N 



2e c _ 



Here, it should be noted that the prefactor x/x s ensures 
that the level density sums up to the correct hole density 
on the stripes. The chemical potential is determined by 
the hole density in the stripes 



= l^ tanh(e - q/ 2T) 



(11) 



For simplicity we assume that the holes act as free elec- 
trons on the lattice with dispersion e q = — 2t[cos(q x ) + 
cos(q y )] - (J,, 

For large critical temperature T c ^t the integral on the 
r.h.s. of the gap equation is simple and we obtain T c — 
Vx/Ax s = Vb(l — 2x/x s )x/Ax s . The pairing strength 
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Figure 2: Critical temperature vs. doping from Eq. lfl3|l for 
Vo/t = 2, 3, 4 from bottom and up, and x a = 1. 



is, however, not strong enough to satisfy that condition. 
For small critical temperature and x s ~ 1 (/x = 0) the 
singularities at e q = allows us to calculate the gap 
integral to leading logarithmic orders 



1 x 



, . - ., [ln 2 (2t/T c ) + 0.41n(2t/T c ) + 1.78] . (12) 
V ir z x s 

Besides the usual BCS log originating from the energy 
denominator another log arises from the logarithmic level 
density due to the van Hove singularity in a 2D lattice at 
half filling. The resulting critical temperature from Eq. 
1H is 



T c = 2iexp 



-y/ir 2 tx s /(xV) - 1.74 + 0.2 



(13) 



The filled stripe condition x s ~ 1 is, however, only valid 
for dopings 1/4^x^,1/2. For lower doping the filling is 
smaller, e.g. x s ~ 0.5 for x < 1/8. The lower density of 
holes generally reduces T c (22| . 

The stripe superfluidity differs from a superfiuid of hole 
pairs. The latter phase would be similar to the molec- 
ular Bose-Einstein condensate (mBEC) discovered in ul- 
tracold atomic 3D traps (without lattices) when a di- 
lute gas of attractive fermionic atoms cross the unitarity 
limit (infinite scattering length) and cross over from a 
BCS to a mBEC (2(|. The critical temperature for the 
superfiuid mBEC is of the order of the Fermi energy in 
3D T c = h 2 {n/2((3/2)) 2 / 3 ir/m ~ 0.218£Vin the mBEC 
and slightly larger in the unitarity limit (35|]. In 2D, how- 
ever, the critical temperature for a BEC vanishes and can 
therefore not explain HTc. 

In Fig. 2 the critical temperatures for dSC of Eq. 
(fl3l) is shown for three effective pairing interactions. By 
construction d-wave superfluidity vanishes near x ~ 1/2, 
where the effective pairing V vanishes, and also as x — > 0, 
where the density of holes disappears. Thus the sim- 
ple stripe pairing model gives a qualitative description 
of HTc. In contrast d-wave superfluidity is largest at 
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Figure 3: Density distributions of Fermi atoms in an optical 
lattice confined by a harmonic trap. U = (e") = 6t and 
x s — 1/2. Full (dashed) curves represents the PS (stripe) 
equation of states (see text). 



low doping for standard uniform Hubbard models with 
an attractive interaction and it extends to large doping. 

As discussed above the semi-phenomenological stripe 
dSC model is just one of many models for HTc dSC. It has 
the interesting consequence that if PS occurs in optical 
lattices there will be no stripes and therefore no dSC 
either. Both will have direct observational consequences 
in optical lattices. 



V. OPTICAL LATTICES IN TRAPS 

We can now address the phases present in optical lat- 
tices based on the equation of states for the 2D Hubbard 
model without Coulomb frustration, and calculate quan- 
tities measured experimentally. 

In the limit of many particles in a shallow confining po- 
tential, the Thomas-Fermi approximation (TF) applies 
because the length scales over which the trapping po- 
tential and density varies are long compared to phase 
boundaries and the lattice spacing [36]. We shall in the 
following assume that the confining potential is on a har- 
monic oscillator form V 2 r 2 as is the case in most exper- 
iments. Within TF the total chemical potential is given 
by the sum of the trap potential and the local chemical 
potential fx(n) = de/dn, where e and n(r) are the local 
energy and number density per site respectively, i.e. 

[i-Tot = n{n) + V 2 r 2 = fi(n = 0) + V 2 R 2 . (14) 

It must be constant over the lattice and can therefore be 
set to its value at the edge or radius R of occupied lattice 
sites, which gives the last equation in (fl4|) . The chemical 
potential for the dilute lattice gas in the 2D Hubbard 
model is /i(n = 0) = — At. 



Measuring the density distribution n(r) gives the 
chemical potential by inverting Eq. (fl4|) . and the equa- 
tion of state can therefore be determined experimentally. 
We can from the above estimates for the equation of 
states also predict the density distribution at least qual- 
itatively. The scale of the chemical potential is rather 
limited since it only varies from /u(0) = — At in the dilute 
limit to a value at half filling between zero for the weakly 
repulsive PM (U < t) and fj,(l) = At for the FM in the 
strongly repulsive limit (U ^> t). Near half filling the 
chemical potential is fx(n = 1) = e± for the stripe model. 
For the PS model the chemical potential is smaller and 
constant between < x < x s , yet on a similar scale 
fi(n = 1) ~ 2t. The average slope of the chemical po- 
tential is therefore (e") = fi{n = 1) — [i(n = 0) ~ 6t. 
Approximating the chemical potential with such a linear 
density dependence [i — (e")n we obtain for n < 1 

nW= (S) (i?2 " r2) - (15) 

The PM and stripe phases be low ha lf filling therefore 
extends up to a radius R c = ^Je" /V 2 . The chemical po- 
tential fi(n) has a gap at half filling which is of order 
U. As result a Mott insulator or AF density plateau ap- 
pears a^. = 1, which extends a radial distance of order 
~ yJU IVi. The phases repeat for densities above half 
filling up to the band insulator at n = 2 such that the 
characteristic "wedding cake layers" appear as shown in 
Fig. 3. The stripe and PS models differ qualitatively at 
densities around half filling where the latter has a dis- 
continuity at n = n s and n — 2 — n s , where we expect 
n s ~ 0.5 as discussed above. The stripe phase is a d-wave 
superfluid between the AF and PM phases whereas no 
such intermediate superfluid appears for the PS model. 
In both cases the AF phase is a MI at half filling n = 1 
only, whereas in the cuprates the AF extends between 
densities 0.95 — 1.2 due to Coulomb frustration or possi- 
bly a finite t' . Detailed measurements of the densities in 
confined optical lattices near half filling could therefore 
in principle reveal the effects of Coulomb frustration. 

The density and radius are related through lfl4|) and 
the corresponding total number of atoms in the confined 
lattice is 

2tt f R 

N = — n(r)rdr . (16) 
a 2 Jo 

The critical numbers for the stripe and AF/MI phases 
to appear in the center of the trap are N = TtR 2 /Aa 2 
and N = nR 2 :/2a 2 respectively. Fillings above half filling 
n = 1 and up to the band insulator BI n — 2 can be 
reached by adding more atoms as shown in Fig. 3. The 
critical fillings and radii depend on the ratio of U/(e"). 
For U » (e") the MI plateau dominates the density dis- 
tribution and the number of particles required to reach 
higher densities and the BI in the center is N — kU /V2C1 2 . 

Other possible scenarios are qualitative different. If 
HTc is described by the 2D Hubbard model without 



7 



Coulomb frustration, then we can expect that the MI 
is replaced by an AF phase around n = 1 surrounded 
by two dSC phases at densities lower (hole doped) and 
higher (particle doped). According to [9fl phases with AF 
and dSC order are mixed when U^&t but coexist when 
U^8t with a first order transition between two densities. 
A critical point is predicted where a mixed AF+dSC 
phase terminate and is replaced by a first order phase 
transition. Here the AFM and dSC phases are separated 
and coexist with a density discontinuity. In []J|, how- 
ever, the pure dSC undergo a first order transition to 
a mixed AF+dSC phase. The densities at which these 
transitions take place correspond to doping x ~ ±0.1. 
If nearest neighbor interactions or next-nearest neighbor 
hopping are included, the particle-hole symmetry is bro- 
ken and the phase diagram becomes asymmetric around 
half filling (n = 1). The models differ dramatically in 
their predictions of the phases in the shells surrounding 
the MI shell. 

Experimentally the density distributions and shell 
structures have been measured for Bose atoms in opti- 
cal lattices by a slicing method in combination with mi- 
crowave transitions that differentiates between singly and 
doubly occupied sites. 0] Very low temperatures and fine 
resolution will be required in order to observe the density 
discontinuities of Fig. (3). 

The momentum distributions of atoms in the confined 
lattice can be found by time of flight experiments. After 
the system has expanded a time t the columb integrated 
densities are measured by light absorption imaging. In 
such time of flight experiments the position is related 
to the momentum in the lattice before free expansion 
as r = kt/m in the far field approximation. As a re- 
sult one measures the momentum distribution (nk) and 
correlations (n-^n^) of atoms in confined lattice. In a 
sudden release (nk) is dominated by the Fourier trans- 
form of the tight binding Wannier functions in the case 
of Fermi atoms whereas a BEC displays characteristic 
Bragg peaks. However, by adiabatically turning off the 
optical lattice potential with Fermi atoms only the crys- 
tal momenta in the lowest Bloch band remains, which for 
a 2D band insulator is a square \k XtV \ < ir, as observed in 
PJ. If the adiabatic release can be made perfect at very 
low temperatures it is in principle possible to measure 
accurate momentum distributions. In a AF+PM mixed 
phase one should observe a two-component momentum 
distribution: a square from the MI and a circular on top 
from the PM. This also requires on-site repulsive interac- 
tion which must be turned off suddenly or corrected for 
in the expansion. 

The density-density correlation function is 



= 5, 



kk'<nk) 



(akak'a k /a k ) 



(17) 



where the first term is the auto-correlations which is sup- 
pressed by a factor 1/N and can therefore be ignored. 
We proceed by expanding the field operator in terms of 
Wannier functions: d(k) = u>(k) J^r e~* kR a(R), over 
the Bravais lattice R. We assume factorization into one- 



particle density matrices (a(R)a(R')) = «,r(5rr/ |39| 
which, however, excludes off-diagonal orders from phase 
coherence in a superfluid. Superfluids will be discussed 
later. In non-superfluid systems we then obtain the cor- 
relation function 



C(k,k') - 



{nkn k >) 
(n k )(n k ') 

N 2 



E^ (k " k ' )R ^ 



R 



, (18) 



with (fik) = Af|u>(k)| 2 . For bosons the minus sign is 
replaced by a positive one and gives the characteristic 
Hanbury-Brown & Twiss bunching observed for photons 
in stellar interferometry and for mesons in relativistic 
heavy ion collisions [40] . For non- interacting Fermi atoms 
the minus sign simply enforces the Pauli principle which 
inhibits two atoms occupying the same momentum state. 

For MI and BI phases of bosons and fermions the oc- 
cupation numbers tir are numbers, and the sum in Eq. 
(HU yields Bragg peaks and dips respectively whenever 
the momentum difference equals the reciprocal lattice 
vectors, i.e. q = k k' = 2ir(n x ,ny), where n XtV are 
integers and in units a = 1. The radius of the occupied 
sites in the optical lattice limits the sum in Eq. (fl8l) and 
results in a finite momentum width of the peaks and dips 
of order the inverse radius. 

Bragg peaks have been observed for bosons in 3D 
3] and 2D [3] lattices, and dips for 3D fermions in 
aj. The Bragg peaks and dips occur in C(k, k') when 
q = Tr(n Xl n y ), where n x ,n y are even integers. In an AF 
phase the periodicity of a given spin is two lattice dis- 
tances and anti-bunching also appears for odd integers 
by an amount [37| C(k, -k) = [(Um/2)/2E k } 2 , where m 
is the AF order parameter and = (Um/2) 2 + e k 
the quasiparticle energy. 

Just as for the AF phase we can in a stripe phase expect 
charge and spin correlations as in low energy magnetic 
neutron scattering, i.e. q = Q c and q = Q s respectively, 
observed as anti-bunching at these wave-numbers. How- 
ever, because the doping x varies in the confined optical 
lattice, the Bragg dips are distributed over the range of 
values for x/x s and are therefore hard to distinguish from 
the background. If, however, the four stripe periodicity 
with x/x s = 1/8 occurs for 1/8 < x < 1/4 as in HTc, 
we can expect novel Bragg dips for charge correlations at 
q = (0, ±7r/2) and q = (±7r/2,0) and for spin correla- 
tions at q = tt(1, 1 ± 1/4) and q = tt(1 ± 1/4, 1). 

Pairing occurs between opposite momenta and leads 
to bunching for k = — k'. Recently, phase coherence and 
s-wave pairing has been observed in optical lattices with 
attractive onsite interactions near the BCS-BEC cross- 
over [3] . As the lattice heights are increased the MI phase 
dominates and phase coherence is gradually lost. For re- 
pulsive on-site interactions the weaker d-wave superfluid- 
ity may appear as discussed above. We can expect atoms 
bunching by the amount C(k, — k) = A k /2i?£, where 

VA k + e£. If dSC is 



the quasiparticle energy is E^ 
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associated with stripes, the bunching due to d-wave su- 
perfluidity should occur in conjunction with the stripe 
anti-bunching. 

VI. SUMMARY 

In conclusion, ultracold fermionic atoms in confined 
optical lattices can reveal details of the equation of state 
and phases of the 2D Hubbard model for a variety of 
interaction and hopping parameters. Specifically two 
plausible models have been investigated with phase sep- 
aration and stripe formation respectively that results in 
qualitatively different density distributions in optical lat- 
tices with and without density discontinuities. However, 
in the cuprates Coulomb frustration can induce stripe 
formation in both cases which makes it hard to distin- 
guish between the two cases. Coulomb frustration could 
as shown in a simple model also be responsible for an 
antiferromagnetic phases extending from densities be- 



tween 0.95 — 1.15 in cuprates whereas in optical lattices 
only ra = 1 MI plateau should form because there is no 
Coulomb frustration. If stripes are required for high tem- 
perature superconductivity, one will not observe super- 
fluidity in optical lattices if phase separation occurs. 

We conclude that the various PM, FM, AF, d-wave 
superfluidity and possibly also stripe phases can be ob- 
served by measurements of densities, momentum distri- 
butions and correlations in expanding ultracold atoms 
from optical lattices. Bunching and anti-bunching cor- 
relations can reveal these phases at characteristic Bragg 
momenta. By varying the filling densities, interaction 
strengths and temperatures, the onset and amount of 
the various phases can be studied in order to determine 
the equation of states. Experimental determination of 
AF, stripe, metal and/or d-wave superfluid phases will 
severely restrict the models, and the density and tem- 
perature dependence on parameters such as U, t, t' will 
provide much understanding of the Hubbard model and 
HTc. 
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